Circulant Hadamard Matrices 

R. Stanley 

An n x n matrix H is a Hadamard matrix if its entries are ±1 and its 
rows are orthogonal. Equivalently, its entries are ±1 and HH l = nl. In 
particular, 

det H = ±n n/2 . (1) 

It is easy to see that if H is an n x n Hadamard matrix then n — 1, n — 2, 
or n = 4m for some integer m. It is conjectured that the converse is true, 
i.e., for every such n there exists antixn Hadamard matrix. 

An n x n matrix A = (bij) is a circulant if it has the form = ai-j 
for some a , ai, . . . , where the subscript i — j is taken modulo n. For 
instance, 

abed 
^ _ d a b c 
c d a b 
b c d a 

is a circulant. Let A = (oj-j) be an n x n circulant, and let £ = e 27n /™, 
a primitive nth root of unity. It is straightforward to compute that for 
< j < n the column vector [1,^, C 2j > • • • > C^ 1-1 ^ 7 ]* i s an eigenvector of A 
with eigenvalue ao + C Ja i + C 2ja 2 + • • • + C^ n_1 ^ a n-i- Hence 

n-l 

det(A) = J](a + C J 'ai + ( 2 >a 2 + ■■■ + C ( "" 1)j a„-i). (2) 

3=0 

Note. The determinant of a circulant matrix is an example of a group 
determinant, where the group is the cyclic group of order n. In 1880 Dedekind 
suggested generalizing the case of circulants (and more generally group de- 
terminants for abelian groups) to arbitrary groups. It was this suggestion 
that led Frobenius to the creation group of representation theory. See [1] and 
the references therein. 
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Note that the matrix 



-1111 

I- 111 

II- 11 

III- 1 

is both a Hadamard matrix and a circulant. 

Conjecture (source?). Let H be an n x n circulant Hadamard matrix. 
Then n = 1 or n = 4. 

The main work on this conjecture is due to Richard Turyn [2]. He showed 
that there does not exist a circulant Hadamard matrix of order 8m, and he 
also excluded certain other orders of the form 4(2m + 1). Turyn's proofs 
use the machinery of algebraic number theory. Here we will give a proof 
for the special case n = 2 k , k > 3, where the algebraic number theory can 
be "dumbed down" to elementary commutative algebra and field theory. It 
would be interesting to find similar proofs for other values of n. 

Theorem 1. There does not exist a circulant Hadamard matrix H of order 
2 k , k>3. 

From now on we assume n = 2 k and £ = e 2m l 2k . Clearly £ is a zero of 
the polynomial Pk{x) = x 2 +1. We will be working in the ring Z[£], the 
smallest subring of C containing Q and (. Write Q(C) for the quotient field 
of Z[£], i.e., the field obtained by adjoining ( to Q. 

Lemma 2. The polynomial Pk(x) is irreducible over Q. 

Proof. If pk{x) is reducible then so is pk(x+l). Recall that by Gauss' lemma, 
an integral polynomial that factors over Q also factors over Z. If p(x), q(x) G 
Z[x], write p(x) = q(x) (mod 2) to mean that the coefficients of p(x) — q(x) 
are even. Now 

p k (x + l) = (x + l) 2 "' 1 + 1 = x 2 "' 1 (mod 2). 

Hence any factorization of Pk{x+ 1) over Z into two factors of degree at least 
one has the form p k (x + 1) = (x r + 2a)(x s + 2b), where r + s = 2 k ~ l and a, b 
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are polynomial of degrees less than r and s, respectively. Hence the constant 
term of pk(x + 1) is divisible by 4, a contradiction. □ 



It follows by elementary field theory that every element m e Z[(] can be 
uniquely written in the form 

u = b + hC + b 2 ( 2 + ■■■ + 6 n/2 _iC n/2 -\ bi e Z. 

The basis for our proof of Theorem 1 is the two different ways to compute 
detif given by equations (1) and (2), yielding the formula 

n-l 

l[(a + Ca, + ( 2 >a 2 + ■■■ + (^a^) = ±n n / 2 = ±2 2 "\ (3) 

3=0 

Thus we have a factorization in Z[£] of 2 2k 1 . Algebraic number theory is 
concerned with factorization of algebraic integers (and ideals) in algebraic 
number fields, so we have a vast amount of machinery available to show that 
no factorization (3) is possible (under the assumption that each aj = ±1). 
Compare Kummer's famous approach toward Fermat's Last Theorem (which 
led to his creation of algebraic number theory), in which he considered the 
equation x n + y n = z n as rir™=i( ;r + T v) = 

We are continuing to assume that H = (a,j-i) is an n x n circulant 
Hadamard matrix. We will denote the eigenvalues of H by 

lj = a + ai ( j + a 2 ( 2j + ■■■ + a n _!C ( "- 1)j . 

Lemma 3. For < j < n — 1 we have 

Thus all the factors appearing on the left-hand side of (3) have absolute value 
y/n. 

First proof (naive). Let Hi denote the ith row of H, and let • denote the 
usual dot product. Then 

Ifij = (a + a 1 C j + --- + a n ^C {n - 1)j )(a + a 1 C j + --- + a n ^C {n - 1)j ) 
= H 1 -H 1 + (H, ■ H 2 )C + (H 2 ■ H 3 )C 2j + --- + (H 1 - H^C,^ . 
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By the Hadamard property we have Hi ■ Hi = n, while Hi ■ H^ = for 
2 < k < n, and the proof follows. □ 

Second proof (algebraic). The matrix -^H is a real orthogonal matrix. 
By linear algebra, all its eigenvalues have absolute value 1. Hence all eigen- 
values 7j of H have absolute value y/n. □ 

Lemma 4. We have 

2 = (1 - C) n/ \ (4) 

where u is a unit in Z[(]. 



Proof. Put x = 1 in 

x n/2 + 1= Y[{x-( 



n-1 



3=0 
j odd 



to get 2 = [1,(1 - C i ). Since 

i-C J = (i-C)(i + C + --- + C^ 1 ), 

it suffices to show that 1 + £ + ••• + is a unit when j is odd. Let 
jj = 1 (modn). Then 



(l + C + ... + C i-i)-i = 



i-C 

1 _ ( C i) 3 



as desired. □ 
Lemma 5. W^e /iawe Z[C]/(1 - C) = F 2 . 

Proof. Let i? = Z[£]/ (1 — £). The integer 2 is not a unit in Z[£], e.g., because 
1/2 is not an algebraic integer. Thus by Lemma 4, 1 — £ is also not a unit. 
Hence R ^ 0. 

For all j we have C J — 1 in R since — 1 = (( — H \~ !)• Hence 

all elements of i? can be written as ordinary integers m. But = 2 in R by 
Lemma 4, so the only elements of i? are and 1. □ 
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Lemma 6. For all < j < n — 1 there is an integer hj > such that 

a + ai (i + a 2 C 3 + ■■■ + a„-iC (n " 1)j = Vj(l - C)\ 
where Vj is a unit in 

Proof. Since 2 is a multiple of 1 — ( by Lemma 4, we have by (3) that 

n-l 

Y[(a + + a 2 C 3 + ■■■ + Gn-iC (n ~ 1)j ) = 

3=0 

in Z[£]/(l — (). Since Z[£]/(l — () is a domain by Lemma 6, some factor 

a + ai( j H h a„_iC^ _1 ^ is divisible by 1 — £. Divide this factor and the 

right-hand side of (4) by 1 — (, and iterate the procedure. We continue to 
divide a factor of the left-hand side and the right-hand side by 1 — ( until 
the right-hand side becomes the unit u. Hence each factor of the original 
product has the form v (1 — () h , where v is a unit. □ 

Corollary 7. Either 70/71 G Z[(] or 71/70 G Z[£]. (In fact, both 70/71 G 
Z[£] and 71/70 G Z[£], as wra// soon become apparent, but we don't need this 
fact here.) 

Proof. By the previous lemma, each 7^ has the form — () h i. If h > hi 
then 70/71 G Z[£]; otherwise 71/70 G Z[^]. □ 

We now need to appeal to a result of Kronecker on elements of Z[£] 
of absolute value one. For completeness we include a proof of this result, 
beginning with a lemma. 

Lemma 8. Let 9 be an algebraic integer such that 6 and all its conjugates 
have absolute value one. Then 6 is a root of unity. 

Proof. Suppose the contrary. Let deg(0) = d, i.e., [Q(0) : Q] = d. Now 9, 9 2 , 
3 , . . . are all distinct and hence infinitely many of them have the property 
that no two are conjugate. Each Qi G Q[0] and so is the root of a monic 
integral polynomial of degree at most d. If 0\, 62, ■ ■ ■ , da are the conjugates 
of 9, then all the conjugates of are among 9{, 9 J 2 , . . . , 9 3 d . Hence each 9^ 
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satisfies the hypothesis that all its conjugates have absolute value 1 (and 9^ 
is an algebraic integer). Thus the rth elementary symmetric function e r in 
$ and its conjugates has at most ( d ) terms, each of absolute value 1, so 
\e r \ < Q- Moreover, e r G Z since 9 3 is an algebraic integer. It follows that 
there are only finitely many possible polynomials that can be the irreducible 
monic polynomials with roots one of the 9 j, s, contradicting the fact that 
there are infinitely many 9^s for which no two are conjugate. □ 

Theorem 9 (Kronecker). Let r be any root of unity and a G Q[r] with 
\a\ = 1. Then a is a root of unity. 

Proof. We use the basic fact from Galois theory that the Galois group of the 
extension field Q(t)/Q is abelian. Let (3 be a conjugate of en, so (3 = w(a) 
for some automorphism w of Q(r). Apply w to the equation aa = 1. Since 
complex conjugation is an automorphism of Q(r) it commutes with w, so we 
obtain @j3 — 1. Hence all the conjugates of a have absolute value one, so a 
is a root of unity by the previous lemma. □ 

We now have all the ingredients to complete the proof of Theorem 1. 
Note that we have yet to use the hypothesis that Oj = ±1. By Lemma 3 we 
have 

|7i/7o| = |7o/7i| = L 

Hence by Corollary 7 and Theorem 9 we have 70 = C~ r 7i f° r some r. Expand 
70 and C" r 7i uniquely as integer linear combinations of 1, (, ( 2 , . . . , C^ _1 : 

7o = a + ai H ha„_i = ±n/2 

C~ r 7i = C~ r ((a - an/2) + (ai - a n/2+ i)C H ) 

= (a r — a„/2+r) + ( a r+l — a n/2+r+l)C + ' ' ' • 

Equating coefficients of (° yields ±n/2 = a r — a n / 2+r . Since each Oj = ±1, 
we must have n < 4, completing the proof. □ 
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